Sensors 2013, 13, 12295-12328; doi:10.3390/sl30912295 



OPEN ACCESS 



sensors 

ISSN 1424-8220 

www.mdpi.com/journal/sensors 

Article 

Currents Induced by Injected Charge in Junction Detectors 

Eugenijus Gaubas *, Tomas Ceponis and Vidas Kalesinskas 

Institute of Applied Research, Vilnius University, Sauletekio av. 9-III, LT- 10223, Vilnius, Lithuania; 
E-Mails: tomas.ceponis@ff.vu.lt (T.C.); vidas.kalesinskas@ff.vu.lt (V.K.) 

* Author to whom correspondence should be addressed; E-Mail: eugenijus.gaubas@ff.vu.lt; 
Tel.: +370-52-366-082; Fax: +370-52-366-079. 

Received: 25 July 2013; in revised form: 21 August 2013 / Accepted: 4 September 2013 / 
Published: 12 September 2013 

Abstract: The problem of drifting charge-induced currents is considered in order to predict 
the pulsed operational characteristics in photo-and particle-detectors with a junction 
controlled active area. The direct analysis of the field changes induced by drifting charge in 
the abrupt junction devices with a plane-parallel geometry of finite area electrodes is 
presented. The problem is solved using the one-dimensional approach. The models of the 
formation of the induced pulsed currents have been analyzed for the regimes of partial and 
full depletion. The obtained solutions for the current density contain expressions of a 
velocity field dependence on the applied voltage, location of the injected surface charge 
domain and carrier capture parameters. The drift component of this current coincides with 
Ramo's expression. It has been illustrated, that the synchronous action of carrier drift, 
trapping, generation and diffusion can lead to a vast variety of possible current pulse 
waveforms. Experimental illustrations of the current pulse variations determined by either 
the rather small or large carrier density within the photo-injected charge domain are 
presented, based on a study of Si detectors. 
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1. Introduction 

The drifting charge-induced current is often a prevailing component in detector signals. An analysis 
of this injected charge drift current is employed for the detection of photons or particles and for 
reconstruction of the electric field distribution over the charge drift area in detectors. However, the 



Sensors 2013, 13 



12296 



problem of the charge-induced currents remains a sophisticated issue in more complicated situations of 
photo-and particle-detectors when devices operate in partial or full depletion regimes with a complex 
field distribution in the presence of carrier capture-generation. Then, the induced current due to the 
charge trapping/de-trapping should be considered. The effects of the acting electric field screening or 
the gain caused by the interplay of the injected carrier charge and the bulk charge of ions in the 
depletion volume should also be involved. As usual, the principles of interpretation of the current 
transients in detectors under injected charge are based on the Shockley-Ramo theorem [1,2]. Ramo's 
theorem was derived based on consideration of Green's theorem and Gauss's law, thus, on the balance 
of electrostatic energy changes and on the analysis of bulk as well as of surface charge induced 
electrostatic fields. Ramo's current is derived at the assumption of the symmetry of an infinite 
electrode and of the spherical invariance of the potential surrounding an elementary charge. Therefore, 
the electrical capacitance of a real system of the finite surface area electrodes and its dependence on an 
applied voltage are out of consideration. Several attempts to generalize a simple Ramo's approach to 
multi-electrode and arbitrary space-charge field distribution systems had been made [3-6], although 
these solutions [3-13] are under debate, as referenced in [3-10]. As pointed out in [7], a simple 
Ramo's relation Equation (19) in [7] between the current within an external circuit and an introduced 
weighting field of geometrical sense is held only in the special case of moving charge in active area of 
a detector with a fixed external voltage applied. Therefore, a simple generalization or direct application 
of Ramo's theorem can be erroneous [7]. Most of the generalization approaches are again based on a 
consideration of electrostatic energy conservation [7-10]. For the space charge systems, like a pn 
junction and pin semiconductor devices, the simple Ramo's relations are not valid, and differential 
weighting potentials [8] are artificially introduced to restore the simple Ramo's type of expressions of 
the convection current. More systematic approaches are employed in [9,10], where space charge and 
different types of charges are included. However, the solutions obtained [9,10] for the estimation of the 
induced current on an arbitrary electrode provide no practical relevance of application with regard to 
semiconductor devices, as pointed out in [10], when the real velocity field of the charge domain is not 
known, and it should be evaluated by means of the complete analysis of field distribution. Thereby, an 
additional kinetic equation should be considered to determine the instantaneous velocity of the charge 
domain. However, applications of Ramo's theorem generally ignore this consideration, rendering the 
expressions obtained not practically applicable. The current pulse shape and duration are directly 
dependent on the temporal variation of the injected charge domain position within the inter-electrode 
space, where the instantaneous velocity and acceleration can be changeable. More complications 
appear within the consideration of a kinetic equation, when the injected charge amount can vary due to 
carrier traps inside the inter-electrode space. 

In this work, the fixed area abrupt junctions with a plane-parallel geometry of electrodes are studied 
by considering the partial and over full-depletion regimes. The vectorial nature of the employed 
quantities of the surface charge domain, of the electric field and of the charge drift velocity is always 
kept in mind. The scalar relations are analyzed along the unite ortho-vector with proper signs ascribed 
to directions. To simplify the understanding of applied models, the single type of induced surface 
charge domain motion is initially considered. This simplification can be sufficient to model current 
transients due to highly absorbed photons or alpha-particles in Si detectors, while the generalization for 
a bipolar charge domain is also discussed. The specific features of the induced charge domain drift 
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currents (ICDC) are revealed within analysis of the simulated ICDC transients and highlighted in 
illustrations of the experimental characteristics, measured on Si pin diodes. 

2. Modelling of the Current Transients Induced by Injected Charge Drift 

2.1. Geometry, Circuit and Electric Field Distribution for Different Regimes 

The same electrical circuit as in Ramo's theorem derivation, is considered: one electrode is 
grounded and the high potential is kept on the other one. These two electrodes are connected to the 
external voltage source in series with a load resistor to register a current transient within external 
circuit, as illustrated in Figure la. The circuit for detection of current transient i(t) routinely includes a 
load resistor and closed input of an oscilloscope type instrument. It is, as usual, assumed that the load 
resistor R L is properly chosen to get the recordable signal, and that a voltage drop on load resistor (iR L ) 
for all the range of detected currents can be ignored, i.e., iRL « U. 

As usual, a quasi-neutral domain of excess carriers is initially generated in particle detectors. It is 
accepted the domains are flat surface vector quantities. These domains injected by light or ionizing 
radiation are characterized by the surface charge density and the direction (vector) of the surface 
normal. Thereby, these domains are directly represented by the electric field of the surface charge. The 
sign (polarity) of the injected charge and the direction of the drift velocity vector are also included. For 
the grounded circuit, the single-side surface charge (and field) is ascribed to the voltage source. The 
drifting domain is also considered as the one- side surface charge vector correlated with drift velocity 
vector direction. The surface charge electric field vectors are initially considered, like the first Poisson 
equation. Then, the scalar equations for an instantaneous field distribution are analyzed. By applying 
an external field source, the injected carriers (by light or ionizing radiation) can be separated into 
oppositely moving surface charge sub-domains q e and qn, which induce charges of opposite sign on 
electrodes. The positive charge on the grounded electrode induced by a drifting charge domain is 
moved by the external source (battery) to the electrode of the high potential and vice versa. The latter 
charge transfer current is actually measured within the external circuit as a signal of either the charge 
domain drift or the charge density change. A sketch of the instantaneous field components is presented 
in Figure 1 . 

External voltage source U produces a positive surface charge a on the high potential electrode, 
which is positioned at the distance d from the grounded electrode. Specific for the junction type 
detectors, the mobile charges within an active volume of a device lead to varied depletion width wo 
dependent on the applied (reverse U r ) voltage on electrode. The grounded electrode is assumed to be 
located at the beginning of the coordinate system (x = 0). The complete neutralization of the depletion 
charge (e.g., eNo + ) in the n-base active layer is obtained through respective depletion (w p +) of the other 
layer (eA^~) of a junction: eNo + wo - eN^Wp + = 0. In the abrupt junction of a pin type diode, it is valid 
w p+ « wo- Therefore, the assumption of the grounded electrode location at x = 0 is valid with 
precision of w p +/wo « 1. The injection of the electron domain, with surface charge density q e , into the 
active volume of a detector at the instantaneous position Xo causes a change of the charge on the high 
potential electrode, in the case of the over full-depletion regime. 
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Figure 1. (a) Sketch of circuit for analysis of ICDC transients. Symbols denote as follows: 
X et h is the instantaneous position of the drifting one-sided surface charge domain of the 
electrons (-q e ) or holes (qh) (with normal vectors n e and n/0, respectively; Ei = E q h + E a + E NDX 
is the electric field caused by the superposition of the hole charge (<?/,) domain, the surface 
(of an area S) charge (+a , with a normal vector n CT ) on the high potential electrode and the 
local surface charge introduced to represent the electric field at x created by a bar (of a 
width (d-xj) of the bulk ion charge (with a normal vector n eNDX )', E2 = E qe + E a + E NDX is 
the electric field caused by the superposition of the electron charge (-q e ) domain, the 
surface charge (+cr) and the local surface charge of the ion bar; E3 = E a + Endx is the 
superposition of field between the electron and hole domains and the surface charge of the 
ion bar; v e ,^ is the instantaneous velocity vector of a drift of the surface charge domain -q e 
or q^ respectively; d is the inter-electrode distance; s and so are the material and vacuum 
dielectric permittivity, respectively; ko is the unit ortho-vector in the uni-directional 
coordinate system; y/h and y/ e are the dimensionless, normalized positions of the drifting 
hole or electron domains, respectively; j CT is the current density; U is the external source of 
voltage; R L is the load resistor. The sketches are illustrated for the electric field 
instantaneous distribution due to the injected negative charge surface domain within a 
partially depleted (b) and fully depleted junction (c), respectively, wo is the steady-state 
depletion width; w q (t) is the charge drift introduced depletion width. 



E 9 ;,=<?j, ££b(-ni,) E =-q e / ee 0 (a e ) 



(U+U FD )/d 
U/d 



, \ 


\ 






n e +a 






-q f 

v. 










' — _ 

E 5 ~ E NDX 






E; = 

Eq e *E 0 + E NDX 



— I— 



—\ — 



(a) 



u 



R, 



U'd 





\ 




"e 






Ei = E*dx 















— t- 



—i H 



(b) 



E^-a'sSai-nJ E =-q e ee 0 (n e ) 



(U+U ro )/d 
U/d 





-q c 











(c) 



Sensors 2013, 13 



12299 



For the partial depletion regime, this leads to changes of the depletion width w q , due to the mobile 
carriers in the electrically neutral region (ENR). As the electrodes are only connected to the external 
measurement circuit, the current induced by the moving charge -q e within the external circuit is 
determined by the surface charge da changes (being a full differential) in time: do/dt. To find these 
changes, the superposition of the induced fields should be considered. 

The junction based detectors contain rather complicated field distribution, and, thereby, need 
specific analysis. The instantaneous field distribution in n-type conductivity region of the p + n abrupt 
junction structure during the monopolar drift of a negative charge in partially and over-depleted base 
region is sketched in Figures lb, c, respectively. 

2.2. Current Transient Induced by Charge Drift in Abrupt Junction 

The analysis can be easily applied to the consideration of the barrier capacitance changes for unit 
area Q, = sso/w q (t). The charge drift dependent depletion width w q (t) should then be introduced, due to 
the injected charge q e . The reverse biased steady-state depletion width wo = [2sso(U + Ubi)/eNDef] m 
serves as an equivalent of the inter-electrode spacing d. Here, Ub\ denotes the built-in potential barrier, 
and Nde/is the effective dopant density. Variations of surface charge a are equivalent to the changes of 
the surface charge at wo, induced by an additional depletion charge bar, as cr~ eN De j(w q (t) - w 0 ). It can 
be proved that consideration of the time dependent changes of the system dynamic capacitance 
Cs q (ascribed to a surface area unite) is equivalent to the analysis of the convection current. In 
semiconductor devices an important role is played by carrier capture and generation processes. Also, 
different regimes of the partial-depletion, of the full depletion and the over-depletion (depending on 
detector width and applied voltage) have their specific features. The mixed regimes of the electrode 
surface charge changes and of the depletion width variations are inherent for the applied voltage values 
close to the values Ufd of the full depletion (FD). For clarity, these different regimes are separately 
analyzed below. 

2.2.1. Current Transients of the Injected Charge Drift in Partially Depleted Detector 

This regime is partially discussed in [14]. Let's consider a regime of the applied reverse voltages 
Um < U < Ufd on the n-type conductivity layer at an assumption that the electron domain is injected 
nearby the metallurgic abrupt junction, and the strength of the electric field there is capable to separate 
the electron-hole pairs, with consequent extraction of holes into p + -type layer. This leads to a 
synchronous change of the depletion widths in the n- and p- type conductivity layers to keep the 
junction system electrically neutral behind the depletion wo, n and w p + width boundaries. To simplify 
the analysis, an assumption of the asymmetric doping of n- and p-layers, i.e., the abrupt p-i-n junction 
is accepted, which enables ones to ignore a voltage drop on p + -layer. It is also assumed that the 
external metallic electrodes are in a rapid dynamic balance with neutral n- and p + - layer material. Thus, 
a rate of the processes within an n-layer region is the slowest one. The latter processes determine the 
current transient caused by a drift of the injected electron domain. 

Using the methodology described in [14], an instantaneous field distribution is obtained by 
integrating the first Poisson equation and by assuming an infinitesimally thin drifting domain of 
surface charge of density q e , as: 
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Here, a vector of the electric field is directed towards the junction, while a surface charge domain of 
electrons can drift towards high potential electrode. To find a depletion width w q and, alternatively, 
Ei(0), due to the injected charge q e domain and the external reverse bias voltage (U r ) source, the 
second Poisson integral should be taken, which leads to expressions: 
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Here, a common depletion boundary condition (E(w q ) = 0) and a proper root of the quadratic 
equation are accepted; /u denotes the carrier mobility. For the reverse biased junction, it is assumed as 
usual [15] that U = U r -Ubu Within a coordinate system at rest, the characteristic time parameters 
TMq,wO and ttof.wo, are defined relative to a steady-state width wo, instead of d. Thereby, the 
characteristic time parameters TM q , w o (time of the Maxwell relaxation of charge q, Mq) and ttof.wo (time 
of flight, TOF) are expressed as: 
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These characteristic times, namely, their equality (T M q, w o = t T of, w o), can be a measure for the validity 
of the electrostatic induction approach. These characteristic times implicate the response time of the 
extended (of distributed charge) electrode (ttof.wo) and of a drifting domain ( rMq.wo)- 

A steady-state depletion width wo is expressed by a well-known formula derived within depletion 



approximation [15] as wo = {2ssoU/eNv) m ■ Consequently, a barrier capacitance is obtained as: 
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where Cm = sso/wo- It can be noticed that w q > wo, and, therefore, the barrier capacitance decreases due 
to the injected charge domain. 
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It can be inferred from Equations (1-6), that the regime of the surface charge domain drift within a 
partially depleted layer of a junction, can only be considered under a few restrictions on relations 
among values of the applied voltage, the doping and the injected charge density. A spatial range for w q 
variations is limited by a geometrical width J of a n-base layer (and consequently by the barrier 
capacitance decrease to its geometrical capacitance value), as: 



2^/ 2a^ 
q pN Def ss Q U \ eN Def (/) 

For electrodes of surface unite area 5=1, this leads to the inequalities, written as: 

85 Q Qe T M,Ndef Id U 

~ — = ~< (1 ) (R) 

eN d ss t 2 w U v J 

Here, TM.Ndef = sso/e/UeNoef is the material dielectric relaxation time within a space charge region. 
Inequality Equation (8) leads to a limitation: 

eN Defd 2 .. U . 

2w Q U FD 

of the injected surface charge which can be moved off by U < U FD . The current density is obtained by 
the analysis of the time dependent variations of a surface charge on electrode due to the extracted 
electrons, as: 

. da TT dC dX e TT dC dw q dX e 

J= = U =U — (]()) 

dt 8X e dt dw q dX e dt {W) 

The rearranged (by these differentiation procedures) expression of a module of the current density 
of the injected charge domain (ICD) drift can be represented as follows: 

I . I 1 1 dX e 1 dX e d Wwn 



2(1 | T TQ F . Wa X 3I2 dt e w 0 z dt e dt (11) 

The obtained scalar form of the current density within a coordinate system at rest (wo) is very 
similar to that of the Ramo's current expression. The main difference is an appearance of a coefficient 
K T , dependent on the dimensionless position y/ = X/wo of a drifting surface charge domain within wo, 
and it is composed of the characteristic times as: 

1 



K = 



2(i+^VT 2 (12) 



The appearance of the coefficient K T is a specific feature of the non-fixed position of the virtual 
electrode (w q ), charge on which is varied by a changed depletion range of ions. Simultaneously, the 
possible drift length is also dependent on the rate of the formation of wo and w q , i.e., on the 
characteristic time TM.Ndef = ssc/e/uNoef = sso/e/un E NR of the stabilization of the transitional /l-thick layer 
(between the depletion and ENR layers) due to extraction of the mobile carriers henr = A^e/from ENR. 
This transitional X layer is related to the Debye screening length [15]. 
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The additional scalar equation (with properly accepted vector direction sign) for a velocity of the 
charge domain drift is now expressed as follows: 



^ = Me E(X e ) = - Me E 2 (X e ) = » e eND «( Wq - X e ) = 1 + - ^) (13) 

The rearranged equation into the dimensionless y/ = X e /wo form can be written as: 
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with the adequate boundary conditions, as t=0 for y/ = y/ 0 and t=td r for y = 1, respectively. This 
Equation (14) is the first-order ordinary non-linear equation, a solution of which is expressed as: 
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Extraction of the y/ (t) function, by integrating Equation (15), might be complicated [16], and it can 
commonly be found by a numerical solution. These solutions y/( t) are only determined for an interval 
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Then, the obtained y/(t) should be inserted into the right hand side of Equations (11) and (12). 
Actually, a direct numerical solution of Equation (14) might be preferable in order to simulate the 
current density transients. The initial component of a rise to the pulse vertex (jicD.FiO) an d the rearward 
relaxation component (jicD,R(0) of a current pulse can also be modelled. Here, it is assumed for 
simplicity that the kinetic equation of motion Equation (14) is only slightly modified due to XMq.wO 
during q e injection. 

For N e carriers located within a domain on its surface area S e , comprising a surface density 
q e = eNe/Se, the injected charge drift is assumed to be a uni-directional process. Therefore, to relate 
more adequately the characteristic times t T of,wo and T Mq , w o, the inequality Equation (9) could be 
rearranged as: 

T TOF,w 0 . 

~ <X (17) 

Mq,w 0 

Thus, for the partially depleted junction, both components d /S e < 1 and (\-U/U F d) < 1 should be 
small, and these conditions ensure that ttof,w(/tmc,,wo < 1. Actually, the correlated drift of the injected 
charge domain (as assumed for the Ramo's regime) can only be implemented at TTOF, W (/ T Mq,wo ~ 1. 
Thereby, the exact Ramo's regime is impossible for the partially depleted junction if K T ^ 1. 
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Generally, variation of an initial component J/caf(0 as a function of time t within the current pulse 
evolution is described for the time interval 0 < t <td r as: 

[1 + ^^/(0-/(0] 
W(0 = ^ ^ (18) 

lT M,Ndef p | r rOF.n 0 *^ f -^j3/2 

Equation (18) describes a component of the pulse with a decreasing current shape within the pulse 
vertex, Figure 2. 

Figure 2. (a) Normalized current density j x TjoF,Jq e transients simulated without 
including of the external circuit impact (red solid curve) and under the impact of external 
circuit with R L = 10 Q. (dash-dotted blue curve). Solid black curve represents the barrier 
capacitance variations during charge drift, (b) The electric field redistribution within the 
partially and the over-depleted diode layer during the drift of the electron domain. 



1.5 




The rearward component {JicdM 0) (Figure 2) of the current pulse is determined by the processes of 
the system capacitance Cb,s q restore to its steady-state value Cb,so = seq/wo- At arrival of the drifting 
charge domain to the wo location, the surface charge field qjsso is equal to the electrostatic induction 
charge field determined by a depletion charge bar eND(w q -wo)/sso. Therefore, the field qjsso is 
completely screened being at wo, as the reverse voltage determined field of positively charged ions 
becomes zero at the point wo- Thus, from the time instant t = tdr, an interplay of carrier diffusion from 
electrically neutral region (ENR) and surface charge field qjsso determines the reduction of the w q to 
the steady-state value wo- This process originates a current (relatively to a domain drift one) 
determined by narrowing of the depletion region, -carriers from the ENR drift into the opposite 
direction relative to the injected domain drift. This current represents the decreasing with time 
(t-tdr) component. 

This current component may be responsible for the appearance of the offset within a current 
transient, inherent for the partially depleted detector. Duration of this process is determined by a 
dielectric relaxation time of the material, namely, TM.Mef- This current flows until the instant tc w o of the 
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barrier capacitance restore to its stationary value sso/wo- It is determined by the barrier capacitance 
charging current. Superposition of these currents leads to a current relaxation component expressed as: 



J ICD,R (0 = 7/czj(^ r )exp[ 



T M,Ndef 



(19) 



applied for the time scale of td r <t< tc w o- Here, C w o = Cb,soS, including the surface area S of electrode. 
As a result, the current density within a pulsed ICD transient is expressed as follows: 



jico(t) = 



L1 + 



(t) = 



1 TOF,w n 



' Mq,w Q 



2r 



for t < t A 



M.Ndef 



[1 + 



' Mq,w 0 



Y(t)f 



(20) 



7/c D ,fiW = 7/ CD (^)exp[ 



(t~t dr ) 

T M,Ndef 



], for t dr <t<t CM 



The offset current relaxation to zero is additionally governed by the parameters of the external 
circuit. The relaxation component of the current pulse is determined by the relaxation processes within 
an RC chain consisting of the system capacitance Cb,s q and a load resistor R L . These elements together 
with the value of the applied external voltage determine the amplitude as well as the shape of the 
current pulse, and duration of the jico.F and }icd,r components. 

2.2.2. Current Transients in Fully Depleted Detector 

The surface charge on a metallic (or heavily doped layer) electrode changes together with the space 
charge bar width due to a moving surface charge domain, when external voltage U is equal to or 
exceeds the full depletion (Ufd) value, U >Ufd- This happens due to a lack of semiconductor material 
(supporting mobile carriers) width (relative to a partially depleted layer case) to withstand the action of 
the electric field due to the injected charge. The equilibrium carriers are extracted into the external 
electrode, if depletion covers completely the active layer width. 

Again, using the methodology described above, a field distribution for the negative drifting charge 
is obtained by taking the first Poisson equation: 
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and the second Poisson integral as: 
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or, alternatively, including the value of Ufd = eNo e jd flsso , it can be re-arranged as: 



Sensors 2013, 13 



12305 



U-U FD = — * + q e —(l-^) (23 ) 

££ 0 £S q a 



The solution for a scalar surface charge density on the high potential electrode is expressed as either: 
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The instantaneous field distribution can then be represented as follows: 



E{xy. 



U~U FD , q e n X ex eN Def d n X e 



E l (X e ) = 7^ + — (l--f) ^Ml--f),for x<X e 

£2( X,) = -^^-ii-^(,-i,f„r *>X, <28) 

■„ d ss 0 d 



It can be deduced from Equations (26 and 27) that similarly to the above considered structures, a 
capacitance of the system is decreased relative to its steady-state value, due to a drifting (X e ( t)) surface 
charge domain. It can be inferred a limitation for the charge density possible to move off: 



^ Ues 0 eN Def d 2 q c eN Def d 2 

e d(X _^e_) W ( i_^) W (29) 
d d 

This q e density also serves for evaluation of the relevant range of the positive capacitance CsNq 
values with q c = UC g = Usso/d. 

A module of the current density, derived from Equations (26 and 27), in the case of the over 
full-depleted (OFD) junction, is again expressed as: 

da 1 dX e 

] °™ = ld =q <-d-d7 (30) 

which formally represents the Ramo's type current component. 

However, the considered situation for a fully-depleted junction is more complicated relative to those 
discussed above. The reason is a degenerated point wfd = d, Figure 2b. From one side, this is caused 
by the synchronous action of the capacitor-specific and the junction-inherent field distribution 
symmetries. As the junction determines a decreasing electric field shape, this field is exactly zero at 
the electrode (wfd = d, the full-depletion condition) and a does not change. Thus, any displacement 
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current component for wfd = d, which is caused by the changes of an electric field behind the drifting 
domain, immediately becomes equal (transferred into) to a conductivity (convection) current, due to 
annihilation of the surface charge domain at X e = d. But the symmetry of a capacitor-inherent field 
implies that the displacement and convection currents, flowing within a capacitor in opposite 
directions, exactly compensate each other to keep the external voltage invariable. The displacement 
current in the over-depleted junction is obtained as: 

1 dX„ 



J displ ^0 



dE 

dt 



[-eN Def d + q e ]-—^,for X e <d/2, (A) 

IdX (3D 

[en 0 d + q e ] e - , for X e >d/2, (B) 

d dt 



This j disp i (Equation (31)) component is different (through an additional eN De jdX/dt type term) from 
that joFD derived by consideration of the changes of a surface charge on the electrode of the capacitor. 
Actually, the current (displacement and drift) components measured within external circuit cannot be 
separated, while the complete current is determined by the changes of charge on the external electrode. 

On the other hand, the displacement current (Equation (31)) eNoetdX/dt and the conduction current 
eriodXJdt components completely compensate each other: the displacement current is caused by a 
space charge bar eN De fif/ e d of ions, and the conduction current enodX/dt (Equation (31)) component, 
which appears due to a seeming extraction of equilibrium electrons n 0 and by producing the surface 
charge enod on electrode, contain the opposite signs in Equations (31). This peculiarity occurs, if the 
space charge bar (eNae/fed) (either behind the injected surface charge domain (when X e = 0) or in front of 
it (when X e = d)) determines an appearance of the displacement current component to exactly 
compensate the conductivity current component enod (to ensure invariance of external voltage U). 
Only for a singular set of boundary conditions U = U FD and X e = 0 as well as X e = d, the displacement 
current ss 0 dE/dt is equal to pE the material's conductivity current, i.e., ss 0 dE/d t= pE, where p 
represents a conductivity of semiconductor material. 

Thus, the complete current within an external circuit is again determined by the displacement 
current due to the injected charge, as obtained in Equation (30). While the space charge bar induced 
displacement currents (within Equation (31)) are employed to exactly compensate the conductivity 
(convection) current component in a fully depleted junction layer. These peculiarities should be kept in 
mind when considering the field of velocities to express definitely the time dependent current changes 
within a current pulse. 

A dimensionless velocity field can be considered by using the accelerating electric field component 
for a geometric width d of the inter-electrode spacing as: 

( . 



\ T Mq T M,NdefJ 



v(X e ) dy/ U 1 eN Def eN Def q 1 1 

— = ^- = // e [— 2+ -if/ + —^y/] = + 

d dt d 2 ss Q s£ Q se 0 d z T0F 2t MNdef 

The coefficients in Equation (32) can easily be rearranged by using the characteristic times 



(32) 
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T TOF 



,1 



£8 



<-M,Ndef 



0 



Ve eN Def 
SS 0 



V e (q e ld) 



(33) 



(34) 



(35) 



and their ratios. It can be noticed in Equation (32), that depending on the relations among values of the 
specific time parameters (Equations (33-35)), and to hold Equation (29) for q e limitation, three type 
solutions are obtained as: 



1 1 

— + - 



^ 0 exp[ — — (— ^ — 1)] + 

T Mq T M,Ndef 

T M,Ndef T Mq 

V (t) = \ Vo + (- — + - )t> 



T T0F ^ T M,Ndef 



{l-exp[-— (— ^ 1)]}, for r Mq >T umf (A) 

_ T Mq T M,Ndef 



1 ID I ^ T M,Ndef 



f° r T Mq ~ T M,Ndef 



(B) 



(36) 



l l 

— + - 



(f 0 exp[ (1- 



l Mq 



c Mq 
T M,Ndef 



)] + 



t TOF ^ T M,Ndef 



{exp[ (1- 



l Mq 

l M,Ndef 



■)]-!}, for r Mq <T Ut 



Ndef 



(C) 



l Mq l M,Ndef 



These equations (Equations (36)) describe the changes of a dimensionless position 0 < y/ < 1 with 
time in the interval of 0 < t < tdr for these three (A, B and C) regimes. The corresponding evaluations 
of a drift time tdr are obtained by using the relevant boundary condition (y/(td r ) = 1) for Equations (36). 
These evaluations of a drift time are expressed as follows: 



l l 

+ 

t tof It 



M ,Ndef 



l Mq 



-ln- 



1 1 

T M,Ndef T Mq 



l Mq 
T M,Ndef 



a-vo) 



i i 

— + 

T TOF ^ T M,Ndef 



-1 



for *Mq > T M ,Ndef ( A ) 



l M,Ndef l Mq 



1 



1 



T TOF ^ T M,Ndef 



for r 



Mq ~ l U,Ndef 



(B) 



T TOF ^ T M,Ndef 



+ 1 



Mq 



) 



1 Mq l M,Ndef 
1 



1 



for *Mq <?M,Ndef ^ 



l M,Ndef 



2t, 



M .Ndef 



1 1 

T Mq T M,Ndef 



(37) 



These different regimes can be realized by varying the applied voltage U (through t T of), the doping 
Noef (TM.Mef) and the injected surface charge density q e (T Mq ). The first regime A (Equations (36 and 
37)) is attributed to a small charge drift. Here, the time dependent variations of the dimensionless 
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position i//(t) are similar to that of the partially depleted junction. These yrft) contain a fast initial 
increase followed by a saturation character of the y/( t) changes when t approaches to tdr- The drift time 
is mainly determined by the injected charge dielectric relaxation time r Mq , modified by a mismatch 
between r Mq and T M ,mef- The second regime B can be associated with a correlated drift of the surface 
charge domain, when y/(t) increases linearly with t, characterized by the invariable tdr, which is directly 
determined by ttof- The third regime C is attributed to the large charge drift, when the injected charge 
is able to locally screen the depletion space charge of ions. Then, ttof = TMq < TM.Ndefi an d the 
correlated (Ramo's type) drift of the injected rather large charge appears. The large injected charge 
determines the increasing drift velocity. Then, neither a drift velocity nor acceleration is constant. This 
leads to an exponentially rising (in time) current density during the charge domain drift time (0 <t <td r ). 
For this regime, the injected charge density is only limited by values q e < qc = C g (l - Ufl/U)U. The 
complete shielding of the external voltage created surface charge a appears if q e > q c = C g U. Only a 
diffusion of the injected carriers is then possible for q e > q c . 

The current density ascribed to different regimes can be modelled by using the relevant expressions 
for \jj(t), taken from Equations (32) and (36), in the drift velocity equation (Equation (32)). The current 
density is then expressed as: 



JoFD (0 " 



1 1 

- + - 



r . 'Mo 'Mo . ,. T Mq ., . t T Mq ... T T0F ^ T M .Ndtf .. , t , T Mq ,.,,„ . ... 

-K— + - — — ) + (l M(Vo ex P[ ( — 1)] + j p-0-expI ( — 1)]}}], for T Mll >T MMef (A) 

T Mq T T0F ^ M.Ndef T M,N<kf 1 Mq T M.Ndef ( 1 L) 1 Mq T M,Ndef 

T M,Ndef T Mq 

9e( + " ). for *Mo =^M,W ; ( fi ) 

T T0F iX MJtki \ JO ) 

1 1 

+ 

i, ,, T Mq T Mq ... 1 Mq .. - t 1 ' Mq ,. T T0F ^ T M,Ndef , . t T Mq ,. .... , ._, 

— [(— + — — ) + G —)\¥a^[ — (1 —)] + —; ^(expl — (1 —)]-!}}], for t tof =T Mq <r Mmf (C) 

T Mq T TOF "M.Ndef T M ,Ndef 1 Mq T M .Ndef ) 1 Mq X M ,Ndef 

T Mq T M,Ndef 

for the respective regimes A, B and C. 

The current density changes within a pulse vertex acquire a relaxation curve shape for the regime A, 
when screening of electrons (drifting charge domain) by ion charge within the depletion width 
prevails. For the correlated screening regime B, a square- wave shape current density pulse appears 
with a flat vertex. While for the correlated (Ramo's type) drift regime C (ttof = TMq), the transient with 
increasing in time current density is inherent. 

The monopolar drift of holes can be expressed using methodology described above for the case of 
electrons drift. The positive charge qn drift is only possible towards the p + -layer. Then, the field for x < X 0 , 
which accelerates q h , is important for the consideration of the induction current. This yields: 

cr = (U-U FD )-±-q h¥h (39) 
a 

The current density, for U = const, again acquires the Ramo's-type expression: 
The drift velocity field is described by a differential equation: 
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d Vh _Mh v _Mh J U 1 eN D d | eN D d 

d 2 £E n ES n 



Vh+ — (n-D 

£Z?n 



Mh T Mq,h 
Me T M,Ndef 



+ 1 



Z Mq,h 1 /4 ^ Mq,h 



T TOF,h 2 // c r Mj 



Ndef 



using the characteristic time parameters defined as: 



1 TOF.h 



' Mq,h 



es, 



M h (<lh /d ) 



Assuming the proper boundary conditions: 

t = 0 for y/ h = y/ 0 
t = t Jr for y/ h = 0 

the solutions of the kinetic equation (Equation (41)) are expressed as follows 

T Mq,h 1 Ml, T Mq,h 1 
T TOF,h ^ Me T M.Ndef 



(41) 

(42) 

(43) 

(44) 
(45) 



/4 T Mq,h t 
<"e T M.Ndef 



■Mq,h 



exp 



///, ^ Mq.h 1 
T M.Ndef 



■Mq,h 



Mh T Mq,h 
Me T M,Ndef 



(46) 



The monopolar drift time is then evaluated as: 



l Mq,h 



-ln- 



2 ^ f r 

T Mq,h 
Me T M,Ndef 



M.Ndef 



+ 1 



Mh T Mq,h j r Mg.ft 1 r Mg.ft ^ 



(47) 



Z^e T M.Ndef T TOF,h ^ Me T M, 



Ndef 



Mh T Mq.h 
Me T M.Ndef 



+ 1 



The current density of the hole drift is expressed (by inserting Equations (41) and (46) into 
Equation (30) (modified for holes drift)) as: 



W)=- 



' Mi/Ji 



I 



Vo ex P 



ft T Mq,h 



+ 1 

Me T M,Ndif , 



exp 



Mh T Mq,h 



+1 

Me T M,Ndif j 



T Mg,h 1 Mh T UqJ> ! ] 



2 ft J" M 



AWef 



ft ^MJW^ 



+ 1 



ft 



+ 1 

^ft ^ M.Ndef , 



l Mq,h 



1 ft T Uq,h 



2 ft 2". 



(48) 



In the case of the positive charge drift within n-base material, holes are always accelerated due to 
acting space charge field. The transient is then observed with the current density increasing with time. 

2.2.3. Impact of Ion Space Charge in Fully Depleted Detector 

The moving charge inside the over depleted space charge layer induces a displacement current 
component, which exactly compensates the conductivity current component, arisen due to a proximate 
contacting of the depleted layer with external electrode (outside layer). As can be inferred for the 
regime B (Equation (38)), characterized by the matched relaxation lifetimes r Mq = TM.Ndef, the space 
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charge eNoef over d accelerates a drift of the injected charge domain by the shortening of the drift time 
to the value t dr = t T ofA^ + (TTOF/2.T M ,Ndef)] < ttof, for y/ 0 =0 (Equation (37)). 

However, for the regime A of the non-correlated relaxation times of the space charge eNoef and of 
the drifting qjd one, existence of the space charge eN De f leads to a reduction of the effective value of 
the drifting charge: 

t 

(49) 



-), for 



T Mq > T M,Ndef ■ 



■M,Ndef 



Then, a surface density of the drifting charge q eye f is instantaneously and locally shielded by the 
space charge of ions, due to the rapid local reaction of the sufficiently large density charge of ions. 
Then, an increase ~ exp[t/T Mq ] of the drift current density (Equation (38)) competes with a seeming 
reduction of the charge density q e ~ exp[-t/ Tm, Ndef], caused by the ion space charge. The space charge 
of ions modifies the current density during q e drift by varying of length of the eNo e fX e bar. 

The large injected charge is able to locally screen the depletion space charge of ions, for the 
regime C. Then, a drift of the injected domain proceeds similarly to that in a capacitor-type device. 

2.2.4. Bipolar Drift of Surface Charge in Junction Structure 



As usual in detectors, a quasi-neutral domain of the excess carriers is initially generated. Then, 
owing to a steady-state applied field, these carriers can be separated into the oppositely moving surface 
charge sub-domains q e and qu- These drifting sub-domains induce charges on the electrodes and 
determine a field in between of them (q e and qh), to hold the initial quasi-neutrality. A sketch of the 
field components is presented in Figure la. 

An instantaneous electric field distribution along the x axis (0 < x < d) for the bipolar drift can be 
represented as: 



E(x) = < 



E l (x) = 


a 


eN D 




ss 0 




E 3 (x) = 


a 


eN D 




ss 0 


ss 0 


E 2 (x) = 


a 


eN D 




ss 0 


S£ 0 



, for x < X h 



(50) 



A field discontinuity at the instantaneous location of surface charge domains is expressed as: 



E = 



E 3 (X h )-E l (X h ) = -^,for X h 
E 2 (X e )-E 3 (X e ) = -^,for X e 



(51) 



Then, the relation between the surface charge +cr on high potential electrode and the external 
voltage U is obtained by taking the second Poisson integral. The solution for a scalar surface charge 
density a can be written as: 



o = (U-U FD )^-q e (l- ¥e )-q h¥h =(U-U FD )C e _ h 
a 



(52) 
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where the expressions are employed for a full depletion voltage Ufd = eNod IIeso and for the 
dimensionless positions y/ e j r = X e> t/d of the sub-domains. However, the full depletion voltage value 
may become variable in the case of carrier capture/emission due to a change of effective dopant 
density N D = N Def . 

It is worth to point out, that in the case of the bipolar drift, the charge cr on the high potential 
electrode becomes dependent on the instantaneous location of both electron and hole separated 
domains of y/ e , y/h). This leads to the appearance of the fields acting on electrons (E 2 ) and holes (£/). 
Theses fields also depend on the instantaneous location of the drift counter-partners as: 



_ U +U FD q e q h eN D d 

E x = — ^- + ^(l-^)-^(l-^ A ) + — 5-^,for x<X h 

d S£ n ££„ ££,, 



U+U FD q e q h eN D d 

d ££„ ££,, ££,, 



(53) 



The induced charge current density, due to a bipolar drift, is expressed as follows: 

. da , IdX, 1 

1 = -* = ^T*- q 'T* ) (54) 



It can be noticed that, owing to v e = -Vh, the scalar current density can be represented by a sum of 
Ramo's-type components: 

d^ + d^ 
H dt dt 

The bipolar drift velocities are correlated during the bipolar drift time n (within time Tb domain) as: 

dy e d ¥h b (56) 
dt dt 

Here, the drift directions are included by accepting the relevant sign for the scalar velocity. There 
exist several situations of the pure bipolar and the mixed drift regimes. These regimes can be separated as: 

= T„, for T b =T dr<e =T iTi 



Wo 


_d-^,)_. 




dWe 


dt 


dt 


¥o 




dWh 




dt 


dt 


W h 


(l-!fo) 


dV„ 


d Ve T 


dt 


dt 



dr,e ~ -dr,h (57) 
for T b =T irM <T ire 

(So) 

for T b = z dr e < T dr h 

(59) 

The regime (Equation (57)) of the synchronous drift of both type carriers within the entire inter- 
electrode gap can only be realized for a single definite point of the charge domain injection y/o. While 
the mixed drift processes appear when the bipolar drift (within current pulse) changes to either the 
monopolar drift of the electron domain after holes reach a p + -layer or it becomes the monopolar drift of 
the hole domain after electrons reach the high potential electrode. These latter situations depend on the 
injection location y/ 0 = y/oh = tyoe within n-base region of the junction and on the mobility of carriers. 
In the general case of the junction structure, the drift process in both layers of the junction should be 
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included. For instance, a drift of holes within n-base region should be extended into p + -layer of the 
abrupt junction to exactly account for the bipolar and monopolar regimes. 

Pure Bipolar Drift 

In the case of the pure bipolar drift regime (Equation (57)), a system of kinetic equations for the n- 
base region and their solutions are expressed as follows: 

dt T b T b (Ol)) 

dy/ e l-y 0 t t 

-p- = — -; v. = v 0 <\ — ) + — (61) 

dt T b T b T b 

These solutions should satisfy the boundary conditions: 

nU=y 0 ; v*U t =o (62) 

(63) 

However, in the more precise approximation, the monopolar drift of holes within p + -region included 
by y/o = y/on + y/op+, should be analyzed. Therefore, in rigorous consideration, the pure bipolar drift can 
be assumed as an idealization. Nevertheless, for the case of T tr ,hp+ « hrjm, the single layer 
approximation can be a relevant approach. 

Then, the inherent time n, of the bipolar drift is obtained by integrating an expression for the drift 
velocity (using Equations (55, 57, 60 and 61)) as: 

dyr h _ n h 



dt 



Ex 



I ^ = i^[4-^^(i-^oa--) + -])-^(i-wi--)]) + ^^a--)]* 

0 d d 2 SE 0 S£ 0 T h T b SE 0 T b SS 0 T b 



WO 



(64) 



¥o T TOF,h 

r b ~ ; ; ; ~ rz ~ ¥o r wF,h 

I 1 ft /l 1 w, \ T TOF,h n , i T TOF,h 

1 + T T0 F,h( )(1-V 0 ) + (l-^o) + - ¥a 

Me T M,Ndef T Mq,e T Mq,h 1 Mq,h 

Inserting these solutions (Equations (60, 61 and 64)) into Equation (55), the current density (in the 
case of the pure bipolar drift) is expressed as: 

] = - = {1 + T — T W F,h( )(l-f 0 ) + (l-ffl) + - ffll* (65) 

T b ¥{) T T0F,h 1 fe T M,Ndef T Mq,e T Mq,h 1 T Mq,h ¥o T TOF,h 

Thus, the pure bipolar drift leads to an invariable current density with pulse duration of t P = y/oT T oF,h> 
Equation (65), provided that carrier capture can be ignored. The dimensionless position of the charge 
injection always is y/ 0 < 1. Therefore, the injected charge drift current pulse is even shorter than a time 
of flight of the counter-partners in bipolar drift, i.e., either t P < ttofm or ttofA^-~ ¥o)- 
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Bipolar Drift During Hole Drift Time 

In the case of the hole drift time x tr ,h is the shortest one among the characteristic times, a bipolar 
(jbB =Ttr,h) drift (Equation (58)) is described by a system of the kinetic equations and their solutions, 
which can be presented as follows: 

dt r M t„ (66) 

^ = 1 ^; ^ e =nexp(-— ) + [l-exp(— ^)] (67) 

These solutions should satisfy the boundary conditions: 

(68) 

V. U=^ 0 ; V. U v ^^;°W, s ) = [l-exp(-l)] + ^ 0 exp(-l) (69) 

Here, ^ e *° serves as the start position for a drift of the electron domain, just during an instant of 
disappearing of the domain of holes at the grounded electrode. The time TbB of the initial bipolar drift 
is obtained by integrating the expression for the drift velocity (using Equations (66-69 and 55)), as: 

d Vh = Mh E 
dt d 1 

? j T ' r n ft, ^ 1 eN D d n 1 u eN D d ft , 9, , , t . , t ,. nl , 

J y h = H^f[--7-~ — — + WoQ )(— — + — )- — <y o exp( ) + (l-exp( ))}}dt 

fO 0 « CI 1 S£ Q T b ££ 0 £S 0 SE „ T trh T trh (70) 

Wo T TOF,h 



T bB - T tr,h ~ 



1 + ^TOF.h 
Me 



lT M,Ndef lT M,Ndef T Mq,e T Mq,e 



Wo T TOF,h 
2 T Mq,h 



Here, the step-like change of field and current density would have obtained for an instant of hole 
arrival to the grounded electrode. To validate the charge, charge momentum, and energy conservation, 
the coordinate transform should be performed, to stitch the solutions obtained in the moving ( TbB time 
domain) coordinate system to that obtained (Equation (30)) for the system of coordinates at rest. This 
transform should include the charge induction on the electrodes, the drift velocity conservation and the 
coordinate relations, which can be represented as: 

a + =a + q 

dt " ( dt dt )U (72) 

V + = \ + (Ve ~ V? e ). *M q ,e < *M ,Ndef 

Y Mq,e 

T M,Ndef 

(73) 



" ¥ + = " i¥e - ¥e e ). > ^M.Ndef 

T Mq,e j 



(74) 



These transforms relate the "new" y/ + coordinate in the system at rest with that y/ (within bipolar 
drift time domain) moving one for the proceeded monopolar drift analysis, after the charge induction 
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procedure is accounted for, and velocities are matched. Then the prolonged monopolar drift velocity of 
electrons (for the case when T Mq ,e < TM.Ndef) is expressed as: 



e,mon y Q,e,mon 



1 r T Ma.e 1 T . 



dt 



Mq,e 1 1 Mq,e 



d T Mq,e T TOF,e 2 ^ M ,Ndef 



*^\n Mq,e 



t M,Ndef 



)] = 



v (),e,mon 



+ Wei + ~ Kl ~ 1 + iYe,mon ~ ¥e )0 " 



0 ' Mq,e 



(75) 



T M,Ndef 



-)] 



Here, and 7? e 2 are the lifetime ratios, which can be dependent on the voltage drop sharing; vo, e ,mon 
is the initial velocity of the monopolar electron drift, which is obtained using relations of velocity 
vectors and their directions for bipolar drift and for the re-calibrated monopolar drift as: 



V e,mon\ [=T ^ — V e,bip V h,bip 



V e,biply/* ,r h 



V h,bip\0,r b 



v Y.bip 



(76) 
(77) 



This gives a coincidence of vo, e ,mon and vm P I (/*" values at the position y/ e *° of an electron domain. 
Expressions for the function y/(t) and monopolar drift time are obtained by integrating (Equation (75)) 
expression with the initial and boundary conditions: 

v,U=vJ°; ^J™„ =1 (78) 

These solutions are given as: 
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(79) 
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T Mq,e 



(1 

The entire duration ?p of the current pulse contains the both phases: 

tp ~ T bB + T e,mon 



(80) 



(81) 



of the bipolar (r^) and the monopolar (r e> mon ) drift of electrons. Inserting these solutions 



(Equations (66-70)) into Equation (55) for the monopolar drift of electrons with y/ e *° as an initial 
position, the current density (in the case of the mixed regime (Equation 58)) is expressed as: 

h=— Wl-exp(-— )] + exp(-— )}, for 0<t<T trh 
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(82) 
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The current transient (for the analyzed case of TM q , e < TM.Ndef) has a decreasing current component 
within the initial phase (during a bipolar drift) and an increasing one within the rearward component of 
the transient due to a drift of the electron domain which screens the space charge. For the case when 
TMq.e = ?M,Ndefi the electrons drift with the constant initial velocity v 0 , e ,mon due to a compensation of the 
drifting charge and the space charge fields. For the case of T Mq ,e > TM,Ndef, the electrons drift with a 
decreasing velocity, as the large space charge of ions (relative to a drifting charge) screens the 
drifting charge. 

Bipolar Drift During Drift Time of Electrons 



In the case the electron drift time r tr<e is the shortest one among the characteristic times, a bipolar 
(Tbc = Ttr.e) drift is followed by a monopolar drift of holes towards p + layer of the junction. After 
performing analogous (as described in previous section) coordinate transformations and solving drift 
velocity equations for the bipolar (during drift time of electrons (Equation 59)) and prolonged 
monopolar drift of holes, the current density is expressed as: 



1h 
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(83) 



with the entire duration t P of the current pulse: 

which consists of the components of the bipolar drift time: 

(1-Wo)t tof , 



The —T, 



and the monopolar drift time: 
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Here, ^ denotes the initial domain position within the monopolar drift of holes; vn,h,mon is the 
initial velocity of the monopolar hole drift which is obtained using the relations of velocity vectors and 
their directions for the bipolar drift and for the re-calibrated monopolar drift as: 
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(87) 

In the case of the proceeded hole monopolar drift within n-base material (after the phase of the 
bipolar drift is finished), holes are always accelerated due to the acting space charge field. The hole 
drift with an increasing velocity determines an inherent shape of the increasing current density within a 
transient, during the monopolar drift phase. 

3. The Impact of Carrier Trapping and Generation 

3.1. The Impact of Carrier Trapping 

The injected charge current can also be changed by carrier trapping and generation. The surface 
charge density dependence on time for the simple traps can be expressed as: 

q e (t) = q e0 CXp(-t/T c ). (88) 

Introducing a trapping dependent dielectric relaxation time as: 



' Mq,tr 



= r M(? exp(f/r c ) 



M e (q M l d)zxp{-tlT c ) 

the time dependent changes of the charge on the high potential electrode are described as: 

eN Dcf (t)d 



(89) 



a(t) = U 



■q e {t)(\-V e (t)) 



(90) 



Here, the time dependent quantities of N De j(t) and q e (t) should be employed. Capture of the injected 
excess carriers, as usual, leads to a synchronous filling of empty donor and acceptor type traps those 
determine the overall charge neutrality and effective doping No e f- Then, the current density within a 
pulse, at assumption that No e j(t)= Afoe/oexpf-z/rt), is expressed as: 
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This equation should be properly matched with the drift kinetic equation: 
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(92) 



The latter equation can only be solved numerically, although the general solution [17] can be 
expressed through complicated integrals [16]. 

A few aspects of the impact of carrier trapping on the injected charge transients for a partially 
depleted junction layer have been mentioned in [14,18]. Evaluation of other parameters (e.g., td r ) of the 
transients determined by the injected charge drift and trapping becomes even more complicated, and it 
can be implemented only by the numerical methods. For prevailing of trapping processes, no 
articulated features of the detector response ascribed to the injected domain drift can be separated, and 
only the relaxation-type shape followed the charge domain injection peak can be observable within a 
current transient. 
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No drift exists (dy/dt = 0) and, consequently, Ramo's current component disappears, if a trapping 
lifetime of the induced charge q e o domain is the shortest one within a set of characteristic relaxation 
parameters. Then, equation (Equation (91)) for current density can be simplified as: 

g eQ exp(-t/r c ) 

J,r (0 = (93) 

T C 

which describes the induction current, - due to the local changes of the injected q e (t) charge density. 

3.2. Generation Current Component 

Carrier trapping, associated with a drifting surface charge domain, may determine the immediate 
(during time significantly shorter than other characteristic time parameters) and local changes of the 
effective charge. A decrease of No e f due to a filling of the charged donor-type traps is equivalent to a 
local charge generation. Thus a simplified approach for evaluation of the generation current can be 
considered. Then, carrier trapping and generation can be analyzed synchronously by rearranging 
Equation (88) as: 

<le (0 = Vet) t eX P (~t 1 T C ) + em Q d eX P (* lT g,ef )] (94) 

Here, mo denotes the initial carrier density on filled traps or the density of the neutralized donor- 
type traps, e is the elementary charge, and r gi ef is an effective generation lifetime. At these simplified 
assumptions, Equation (94) can further be rearranged as: 

2 e (0 = <leo exp(-f /r c )[l + g m0 exp(t/r gCef )] = q e0 exp(-t/r c )f(t) (95) 
with the additional designations as: 

e m o d (96) 
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f(t) = [l+g m0 exp(t/T gCef )] . ( 98 ) 



gC ' ef >c W ' (97) 



This simplified approach enables one to include into consideration the local charge generation. 
Unfortunately, the solutions can be obtained only by numerical analysis. 

The discussed above simplified analysis of the components of carrier drift, trapping and thermal 
release enables one to make the rough estimations of the impact of different components. However, the 
rigorous consideration of processes should be based on the causality principle. The current changes can 
only appear during or after injection of q e . Therefore all the relations for the electric fields and charges 
caused by q e , obtained within the electrostatic approach, and containing the time dependent q e (t) 
components should appear as the convolution integrals, for instance, as: 

j(f) = — = — l^^(l-y(t-&))d® + ^ 1 f? e (®)-^ — d® ; Kyy) 

dt t l0 t c 2t c t l0 dt 



Sensors 2013, 13 



12318 



^_J_|exp[-a-0)/, c ]^_ QVQ _J_|exp[- (? - Q )/, c ] JQ _J_ = o 

0 T M,Ndef 

(t-Q) 

T TOF 

This specification leads to the integral-differential equations. These equations can only be analyzed 
numerically. Then, the above presented simplified models can be employed for the initial and 
qualitative prediction of the numerical solutions. 

4. The Role of Carrier Diffusion 

The drift velocity is varied through the electrostatic interaction of injected charge q e and surface 
charge a on electrode, due to external voltage source. The initial zero drift velocity actually appears 
when the detector signals are caused by the secondary electron-hole pairs generated by the energetic 
elementary particles (high energy photons, hadrons, etc.). Then, a neutral domain with an equal density 
of electrons and holes in pairs is locally generated. The external field is able to separate and move 
these counter-partners towards opposite directions if the density of these carriers is less than q c , i.e., 
for q e < q c = C g U. 

In the partially depleted junction layer, for U < Ufd, the current (ascribed to the injected charge 
drift) varies due to the temporal changes in w q (t), and it really contains a pair separation X e .h(t) = X e -Xf, 
length. The charge separation process induces the change of the depletion width w q (an increase, 
relatively to its steady-state value wo, n & P +) in both layers of the junction, i.e., w q ,„ and w q , p+ . The 
extracted excess holes are located at p + -side producing the same value of the surface field. Thus, the 
overall charge balance w p + N~a, p + = WnN*D, n together with q hiP + X 0iP + = q e , n Xo, n is maintained in a 
diode for the moderate density of the injected charge domain. In the heavily doped layer, the 
characteristic times, e.g., dielectric relaxation times, are significantly shorter than those in the resistive 
layer. Therefore, current transient pulse duration is mainly determined by the longer processes within 
more resistive layer. This motivates an approach of a separate consideration of electron drift within a 
base region of the reverse biased diode detector, where a charge separation process is assumed to be 
sufficiently short, and a drift starts after extraction of counter-partners (separation of pairs) process is finished. 

In the general case of the local injection of excess carrier pairs, separation of counter-partners 
depends on their densities. The external source induced charge <j on electrode can be completely 
shielded by the large injected charge q e within a Debye screening length during the dielectric 
relaxation on metallic electrode, which is extremely short. The space charge of ions is also screened 
during TM q , which is then the shortest one among the characteristic times, in the case of q e » eNo e jd 
and (q e /sos)d > U. The external source is able to react by changing cr, till the system capacitance is 
Cs q > 0. However, the large injected charge q e nearby the grounded electrode (X 0 = 0) is able to create 
an internal field and a voltage drop (q e /sos)d, which reduces the dynamic capacitance of the system to 
Cs q = 0. Then, the external voltage source is completely blocked in supporting of a. As a result, no 
separation of the electron-hole cloud (into domains of electrons and holes) appears. Thereby, 
relaxation of the injected quasi-neutral domain happens completely by carrier diffusion process. 

The reason is the excess carrier diffusion and appearance of the diffusion induced inner field [19] as: 

E D = 1 (D h Vp-D e Vn) (101) 

PMh +n Me 
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This field is proportional to the excess carrier density gradients. For the locally generated domain of 
the nearly infinitesimal width (for instance in tracking of hadron path) within significantly wider inter- 
electrode gap of detector, diffusion due to a sharp gradient (which is also proportional to the carrier 
density) induces an inner electric field which balances a further widening of the domain. Strength of 
this field can be sufficient to compensate partially or fully the applied external field (surface charge cr) 
in the case of low applied voltage and the rather high densities of the injected carriers. Then, carrier 
diffusion and drift leads to the outspread of the domain. Such a process is characterized by the 
ambipolar diffusion coefficient D a . In thin detectors with small applied dc voltage, a width of the 
outspread domain can approach the geometrical dimensions of the inter-electrode layer, even when it is 
formed by the strongly absorbed radiation. The same situation can be realized by the photo-injection of 
excess carrier into a rather thick detector using the homogenously absorbed light. Then, the injected 
domain sweeps the inter-electrode spacing. The external electric field acts as the accelerating factor for 
the surface recombination sd/D (of velocity s, related to sd/D = d/L D by Debye screening length L D ). 
This problem is very similar to the excess carrier ambipolar diffusion moderated carrier recombination 
on surfaces. Then, action of the external field and carrier drift can be included into the properly modified 
surface recombination velocity s. Solution of this problem is well-known [20-23], and it leads to a 
transient of conductivity current due to the carriers arrived to the electrode and extracted to its surface. 
Time variations of the excess carrier density in this domain (averaged over the geometrical thickness d 
of inter-electrode space) is expressed through a sum of the space mode r\ m components as: 

CO 

n ex (0 = n ex (t = 0)£ A m cxpi-rjlDj) (102 ) 

m=l 

The space frequencies ( rj) of these decay modes are described by the solution of the transcendental 
equations of type: 

, s, r Dri s ? . 
C ^ = WT~W (103) 

Here, s } = L D j/r T oF,LD,h and S2 = L D Jr T0 F,LD,e denote the surface recombination velocities ascribed 
to the electric field caused extraction of carriers towards the front and rear (surfaces) electrodes (with 
relevant Debye lengths L D > L D:e ), respectively. These (Equation (102)) solutions with roots found 
from Equation (103) predict a two-componential, the relaxation type current transient. Such a transient 
contains the initial non-exponential relaxation component. The asymptotic decay component is 
characterized by the time parameter td = l/r/j D of the main decay mode, representing the effective 
time of the domain dissipation, — due to diffusion over the entire inter-electrode spacing. 

5. Current Transient Changes Determined by a Signal Recording Circuit 

A signal registration circuit (namely, load resistor) inevitably transforms the current transient shape. 
This appears due to the voltage sharing and the consequent change of a voltage drop on detector 
depending on current value within the circuit. In more general case, the transients are described by the 
solutions of the differential equation with variable coefficients, derived as: 
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This leads to a differential equation: 

di d\nSC s Jt) 



U d\nSC s Jt) 
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dt dt 



R L SC s Jt) 
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which should be solved by using the initial conditions: 

i 0 (t = 0) = 0 
SC Sq (t = 0,q e =0) = C, 



o 



(107) 



(108) 



for the ascending component of a transient, and: 



hAt = t dr ) = j(t = t dr )S = i(t dr ) 
SC Sqir (t = t dr ,ys e =\) = C 0 



(109) 



(HO) 



for the relaxation stage of the transient, respectively. 

The changes of a system dynamic capacitance determine the initial delay and the final stage 
(relaxation) components within the simulated transient. These components are inevitable within the 
charge drift current transients, recorded in experiments. Also, these components should be included 
into the evaluation of the charge collection efficiency. Depending on the geometrical capacitance (C g ) 
and load resistance (Rl) values, the current pulses are significantly modified. 

6. Discussion 

The simplified models [23,24], based on Ramo's expression for the drift current, are attractive as 
they provide a simple analytical description of the detector signals. However, the analytical 
expressions can only be obtained for the simplest approximations. The analytical form of the correlated 
drift (Ramo's-type) current for the junction type detectors is only applicable for a primary estimation 
of a transient shape. Different regimes in the formation of the pulsed response of detectors can appear 
in a real measurement. The time-dependent variations of the current transients may be determined by 
the injected charge dissipation through the domain drift, dielectric relaxation (due to media 
polarization effects), through carrier capture and thermal release processes in the traps containing 
material, via ambipolar diffusion processes. Several specific aspects of these phenomena have been 
discussed above. 

6.1. Limitations of Models 

Adaptability of the simplified models presented above is additionally limited by several factors. A 
principal limitation leads to the threshold values of the acquired drift velocity that should be 
significantly less than those of the electric field (light) propagation velocity in the material under 
consideration, to ensure the validity of the electrostatic approach. This condition excludes the 
possibility to detect the primary charged particles (moving with relativistic velocities) within the 
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inter-electrode spacing. Thus, the secondary particle (the electron-hole pairs with a zero initial drift 
velocity at the injection point) induced currents should be calibrated to the primary particle impact. 
The specific feature of the prevailing drift current caused by the monopolar charge domain is the 
increase of the drift current with time within the vertex of a current pulse. To separate the neutral 
domain (locally generated) into the drifting charge sub-domains, the sensitivity threshold for the 
applied voltage appears. This limitation leads to a condition of the elevated values of bias voltage, at 
least U > Ubi for the partially depleted semiconductor detectors. An enhancement of the external 
voltage may lead to the carrier velocity saturation, which complicates the analysis of the drift velocity 
field: Vdr(X e ). Values of the highest external voltages are also restricted by the necessity to exclude 
repeated and non-linear drift processes, — such as the photo-electric gain moderated by carrier trapping 
and the avalanche processes of the impact ionization or Pool-Frenkel effect. 

6.2. Limitations in the Evaluation Precision of the Depletion Layer Boundary 

In the analysis of the junction type detectors, the parabolic approach has been employed, which 
relates the applied voltage and the width of the depleted region, and it is routinely exploited in device 
physics [15,25]. This approach enables one to simplify the expressions for the electric field 
description, as discussed in [15]. But this approach limits the precision in evaluation of the 
characteristic widths of wo, w q and wfd- In the steady-state case, the Debye screening length L D can be 
a measure for the evaluation of the precision in the estimation of wo- Thus, the transition layer [15] 
width, expressed as: 

2S£ Q 1/2 1/2 

A= W Ux ' cm) 

depends on a voltage drop Ux, ascribed to this transition layer, and on the effective doping of a 
material. Thereby, for the large resistivity material, a relative inaccuracy (A/wo) of the determination 
of the effective widths w 0 , q ,FD ± A. can be unacceptable in the range of small reverse bias voltages. 
Traps and their filling processes can also be a reason for the instantaneous and local changes of the 
effective doping. Thereby, the resolution limit in time scale is expressed through the relaxation times as: 

_ T TOF,U ' x 

T M,Ndef - ~ (112) 

with ttof, ux = A 2 /jU e Uz. The latter condition is determined by a necessity to stabilize the geometrical 
boundary by the balance of the local fields of carrier diffusion and drift. 

The transitional layers are actually inherent to the boundaries between the metallic electrodes and 
dielectric or external heavy doped layers of junctions. Owing to a short dielectric relaxation in the 
heavy doped layers, the semiconductor junction is preferential relative to a dielectric in between 
of electrodes. 

6.3. The Impact of the Injected Charge Density 

In the Ramo's derivation of the charge drift current, it was clearly proved the reciprocity principle: 
the reversibility and equality of the mutual action and reaction of the charged electrode and drifting 
charge. This is based on the conservation of the charge (cr induced electro-statically on the electrode 
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plane by a drifting injected charge q = So), charge momentum qdy//dt (qvdr) and the electrostatic 
energy (q0 = o U). Here, 0 is the surface equi-potential surrounding the drifting charge. The main 
equation for the current density (for instance, considering the motion of the electron domain) can be 
directly obtained by using the electrostatic energy balance Sq0 = ~S(oU). Here, S means a change in 
the electrostatic energy due to a variation of the surface charge (Sa) on the electrode, which should be 
balanced by a change in the energy of the moving charge q e . In the latter balance, q e is assumed to be 
invariable, and these energy changes are ascribed to the changes of the potential S<P(X e ), — during 
charge drift. The temporal changes of the surface charge on the electrode gives current density variations 
dependent on time (for a fixed external voltage), and this current density is generally expressed as 
j(t) = do/dt = -(q/U)(d0/dX e )(dX/dt). Accepting the general electrostatic relation 
E= -gmd<P = n e q e /sso = E q and assuming that instantaneous E q (y/ e ) = q e y/ e leeo, for its scalar 
representation, the expression for the current density is rearranged as j(t) = do/dt = 
q e (d 2 /U)(q e /dsso)[d(X/d)/dt] = q e ( TTOF,/TM,q)[d(X/d)/dt]. Thereby, the just derived current density (on 
the basis of electrostatic energy conservation in the case of our consideration) is consistent with 
Ramo's derivation (also made on the basis of energy conservation) if ttofJtm^ = 1. On the other hand, 
the equality of TroF.e and Xm a is consistent with electrostatic induction approach. Using the scalar values of 
the field within the inter-electrode space as E a = Uld and divE q = VE q ( y/ q ) = q e /ssod, for the injected 
charge field E q (over a geometrical width d), and a balance of electric fields and 
E a = Uld = divE q = V-E q f y/ q ) = q e /ssod, one gets a weighting field We = divE</E a = d~ l . This result 
validates the equal action and re-action of the surface (o and q) charges and leads to the equality of the 
response times tM q , e = TroF.e- Actually, the equality of the response times TMq, e = TroF.e is ensured due to 
correlated changes of the acting voltage which is varied with y/(t). To find the drift velocity field 
Vdr(X e ) = dXJdt, the problem should be solved by consideration of the fields and charges in details. 

As it has been demonstrated above, the drift velocity Vdr is a function of the instantaneous charge 
domain position and the characteristic times: TroF.e, Tm, 9 and tM.Ndef. It can be proved that the current 
density j(t) = 9s£ 0 juU 2 /8d 3 obeys Mott-Gurney's law [26,27], for v dr (X e = J/2) = ju[(U/d) + (q/2s£ 0 )]. 
Therefore, all the applied voltage U drops within the gap between the high potential electrode and the 
q e domain, during the initial phase of a drift process. The drifting domain additionally acts as a voltage 
sharing element (divider) with the parabolic-like characteristic of U2 = (l-yf)U. As a consequence, 
the drifting charge coordinate y/(t) ~ exp(t/TM q ) increases exponentially with the drift time t, leading to 
a variable drift velocity Vd r ~ dy4t)/dt) ~ (d/TMq)exp(t/TMq) and the acceleration a(t) ~ (d/tM q )exp(t/TMq), 
to hold the processes (of drift, voltage sharing, and the induced amount of charge) correlated in time. 
In the case of the pure bipolar domain, the correlated drift of q e and q% sub-domains is equivalent to a 
widening of the quasi-neutral e-h domain. 

In real detectors, the prevailing regime is the detection and collection of a small drifting charge, 
where XMqJtTOF.e » 1. Unfortunately, this regime can only be approximately considered within the 
one-dimensional approach. The reason is a slow dielectric relaxation (xMqJtTOF.e > 1) of a small 
drifting charge q e . Due to the small q e , the charge domain surface becomes corrugated under the action 
of the electrode charge (at boundaries of the detector electrode of finite surface area S) and the charge 
density gradients within the surface plane of a drifting domain. To stabilize the gradients (or the 
oblique action of o surface segments), the drifting charge should vary its position in all three spatial 
dimensions. Thus, the lateral fields should be taken into account, — the charge of dangling bonds on 
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perpendicular boundaries acting as the surface recombination sinks can be a reason for such fields. On 
the basis of the Lagrange variational principle, it can be understood that charge movements within both 
the electrode and the domain planes should be correlated, to react most rapidly on each others changes. 
Then, energy conservation can only be considered by an analysis of the three dimensional charge drift 
and diffusion problem. This leads to the appearance of charges and their neutralization currents on the 
perpendicular (to an inter-electrode drift direction) boundary planes of the base region of the junction 
type detector. The small drifting charge q e is able to terminate the electric field of the electrode's equal 
to its amount. Therefore, a small drifting charge acts as a linear voltage divider within the 
inter-electrode gap, and, consequently, the I-V characteristic obeys Ohm's law. Then, due to the charge 
drift, the acting voltage is different from that applied on electrode U, and it changes under the variation 
of the charge position. The small charge drift current density is less (in comparison with the regime of 
the large charge drift) due to the less effective voltage U e f, acting on the injected charge, q e . This leads 
to the increased TroF.e ( t). The approximation models for U e f should be applied to cover the entire range 
of the injected charge values. Namely, for an electron domain the ratio R e i = TMq,J^TOF,e and 
R e 2 =TMq,</TM,Ndef should be replaced by R ef>el = (T M q,e/rTOF,e)(U e fJ(U + Ufd)) and 
R e f,e2 = (^Mq,e/^M,Ndef)(U e f,</(U + U FD )) (which modifies the voltage) using the approximation 

1/2 ^0 

U e f, e = l.225q e d( y/ e ) /sso. Here, y/ e is the dimensionless location of a drifting domain at the end of 
a bipolar drift. While, for a hole domain drift, these are R e f,hi = {^Mq,hl^TOF,h){U e f,h/(U - Ufd)) and 
R e fM = (TMq,h/ T M,Ndef)(U e f,i/(U - U fd)) (which modifies the voltage) using the approximation 
Uc.ef.h = l.995qhd(y/h*°) lB /££o- These approximations can be understood by an equal redistribution 
between degrees of freedom for the three-dimensional motion, if R e> h ^ 1. The applicability of the U e f 
approximation models has been verified by their relevance to stitch together the one-dimensional 
solutions of the bipolar and the monopolar drift, thereby matching the synchronous conservation of the 
charge, charge momentum and current density continuity. This enables one to get continuous current 
density variations within the simulated vertex of the charge drift current pulse. 

6.4. Correlation with Experimental Results 

As mentioned above, a vast variety of possible pulsed current transients, composed of drift, 
diffusion and displacement current components exists depending on the detector design and different 
external factors, such as the injected charge quantity, applied voltages, presence of traps, etc. 

The simulated specific transient shapes associated with different regimes of the injected charge drift 
are illustrated in Figures 3 and 4. These simulations have been performed by using the above presented 
models and keeping nearly the same charge drift conditions, while varying external voltage to 
implement the partial or full depletion regimes. The described analytical solutions enable ones to get 
the continuous curves of the drifting charge velocity and of current density as a function of a domain 
position and of time. These illustrations demonstrate that current density transients, containing a rather 
flat pulse vertex, can be found in experiments implemented for the small injected charge density by 
using a small load resistance. However, only the largest current components can be resolved due to a 
weak signal. The double peak containing current transient shape should be observed as usual in the 
case of the relative large charge drift, Ramo's type, regime. Thereby, depending on the applied load 
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resistor and voltage, a vast variety of current transient shapes can be obtained in modelling of 
detector responses. 

Variations of current transients due to the injected charge drift, observed in experiments, are 
illustrated in Figures 5. These variations are demonstrated as a function of the injected charge density, 
in Figure 5a. These transients were recorded on the non-irradiated CERN standard pad-detectors, made 
of pure Si of 2 kQ cm resistivity p + nn + structures. The reverse bias voltage was kept fixed with the 
rather moderate values U < Ufd- The surface domain was injected by strongly absorbed light laser 
pulse of 400 ps duration. Density of the injected domain of electrons, initially located nearby a 
junction boundary, was varied by changing intensity of excitation laser beam. The laser pulse was 
sufficiently shorter than RC = 1 ns of the measurement system. The pulsed current was detected on 50 
Q load resistor and registered by a 1 GHz band real time digital oscilloscope. 

Figure 3. Transients of the normalized current density j x TTOF,e/q e simulated for the large 
(a) and small (b) monopolar charge drift, simulated using different values of the carrier 
capture lifetime re- Solid curves are calculated using the analytical expressions, while the 
intermittent curves are obtained including the RC of a signal recording circuit. 
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Figure 4. Transients of the normalized current density jx Zrovjqe simulated for the bipolar 
charge drift during hole (a) and electron (b) drift times for various regimes, dependent on 
the densities of the drifting charge and of the space charge of ions. Solid curves are 
calculated using the derived above analytical expressions, while the dotted curves are 
obtained including the RC of a signal recording circuit. 
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The transients, — characteristic to the pulse durations controlled by the ambipolar diffusion lifetime, 
are illustrated in Figure 5a using a logarithmic time scale. Here, the pulse duration is varied in the time 
scale from a few of nanosecond, — that is inherent for the electrons drift time in the base region, to a 
few microsecond of the diffusion time. It can be easily noticed that current increases with time within a 
vertex of a pulse for the smallest injected charge densities, as predicted by models described. The 
current pulse has nearly exponential relaxation component, after domain reaches the rear electrode. 
Duration of the pulse vertex, measured between the initial (which is on the left in this scale) and 
rearward kinks within current transient is well correlated with electron drift time in d = 300 um thick 
Si layer, using ju e = 1,220 cm /Vs value of the electron mobility. The enhancement of the injected 
charge density, proportional to the n ex o injected carrier concentration, leads to the increasing delay (of 
the rear kink in current transient) and to an increase of the current (proportional to q e ), approaching to 
the ambipolar diffusion lifetime td = d 1% D a , ascribed to the main decay mode. The extracted value of 
the D a = 15 cm /s using the measured td time is in good agreement with parameters ascribed to a 
rather pure Si material. 

Figure 5. (a) Variation of current transients, registered in Si pad-detector biased with a 
fixed reverse voltage of value U < Ufd, as a function of the injected charge density, 
(b) The current transients, ascribed to the charge domain injected by a laser pulse, as a 
function of 25 MeV neutron irradiation fluence, evaluated from the exposure time of the in 
situ experiment. 



Variation of current transients measured at extremely small excitation densities (close to those 
possible to detect at a threshold sensitivity of the measurement system equipped with proper current 
amplifier) is illustrated in Figure 5b. These transients were recorded during 25 MeV neutron irradiation 
when density of radiation induced traps varied with neutrons fluence. The electrical (U r = 150V) and 
optical (n ex o) parameters were carefully controlled to be fixed within measurements. The sample was 
kept in air just behind the neutron beam cone, while other experiment details are published in [14,28]. 
Evolution of the current transients is illustrated in Figure 5b, where currents had been controlled 
starting from that registered in the non-irradiated diode up to exposure duration for which the collected 
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irradiation fluence reaches value of > 10 n/cm . The transient waveform inherent to the drift 
dominated current (curve 1 in Figure 5b) is observed in the non-irradiated diode, which is coincident 
with modelled transient shape. The radiation introduced traps determine a rapid reduction of carrier 
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lifetime and an enhancement of carrier capture and generation current components. The transient of the 
carrier capture dominated current contains the single peak, and the relaxation-type curve occurs 
(curve 4 in Figure 5b). In the range of the intermediate exposure time (curves 2 and 3 in Figure 5b), the 
current pulse duration sustains values of the electron drift time within a diode base. The double peak 
and rising pulse vertex shapes alternate during increase of fluence. Applying of the smallest possible 
densities of the surface charge was sufficient to hide the inhomogeneities of the photo -generated 
domain. This evolution of current transients can be explained by the competition of carrier 
capture/generation and drift currents. For the largest irradiation fluence, the carrier trapping process 
becomes dominant, while a drift current component is completely competed, and the relaxation-type 
current pulse is only observed. 

7. Summary 

The models of the formation of the injected charge pulsed currents have been developed concerning 
the junction-type detectors. The partial and full depletion regimes have been analyzed. It has been 
shown, that, in junction detector, the drift time for the rather small density of the injected charge is 
shortened relatively to that of the capacitor-like detectors when a proper frame of reference (for 
comparison) is accepted and the characteristic relaxation times are matched. The description of the 
current pulse shape for the large injected charge drift in a finite area detector is coincident with that 
derived for the correlated drift (Ramo's-type) expressions. However, the induced currents obtained for 
the regimes of the small injected charge and of partial depletion lead to deviations from the Ramo's 
expressions. The analysis of the drift velocity field revealed the current increase within a vertex of the 
current pulse, for the monopolar drift regime. It has been shown, that presence of carrier traps 
considerably modifies the shape of the current transients. For the extremely large density of the 
injected charge q > C g U, the ambipolar diffusion of the injected carriers may become dominant in 
formation of the injected charge current pulse. It has been illustrated, that synchronous action of carrier 
drift, trapping, generation and diffusion lead to a vast variety of possible current pulse waveforms. 
Experimental illustrations of the current waveform variations obtained for both the rather small and the 
large charge density of the photo -injected domains are presented, based on study of Si detectors. 
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